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ABSTRACT 

We present new more sensitive high-resolution radio observations of a compact broad absorption 
line (BAL) quasar, 1045+352, made with the EVN+MERLIN at 5 GHz. They allowed us to trace the 
connection between the arcsecond structure and the radio core of the quasar. The radio morphology 
of 1045+352 is dominated by a knotty jet showing several bends. We discuss possible scenarios that 
could explain such a complex morphology: galaxy merger, accretion disk instability, precession of 
the jet and jet-cloud interactions. It is possible that we are witnessing an ongoing jet precession in 
this source due to internal instabilities within the jet flow, however, a dense environment detected in 
the submillimeter band and an outflowing material suggested by the X-ray absorption could strongly 
interact with the jet. It is difficult to establish the orientation between the jet axis and the observer in 
1045+352 because of the complex structure. Nevertheless taking into account the most recent inner 
radio structure we conclude that the radio jet is oriented close to the line of sight which can mean that 
the opening angle of the accretion disk wind can be large in this source. We also suggest that there 
is no direct correlation between the jet-observer orientation and the possibility of observing BALs. 
Subject headings: 



1. INTRODUCTION 

Broad absorption lines (BALs) - high ionization reso- 
nant lines (CIV 1549A) and low ionization lines (Mgll 
2800A), are seen in a small number (15%) o f both the 
radio-quiet and rad io-loud quasar populations (jDai et al.1 
20081 iKnigge et all 2008; G ibson et al l 120091). according 
to the traditional BAL QSO definition of lWevmann et all 
(1991). They are probably caused by the outflow of 
gas with high velocities and are part of the accretion 
process. Th e presence of BALs is a geometrical effect 
(lElvisll2000D and/or is connected with the quasa r evolu- 
tion (jBecker et~al1[2000): IGregg et al. 20 0(1 EfflOfih. 
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Theoretical models (jElvisI 120001; Murray et all 119951) 

suggest that BALs are seen at high inclination angles, 
which means that the outflows from accretion disks are 
present near the equatorial plane. However, some re- 
cent numerical work indicates that it is also plausible 
to launch bip olar outflows from the inner regions of a 
thin disk (e.g. iGhosh fc Punslvl[2007t iProga fc Kallmanl 
2004). There is growing observational evide nce indicat- 
ing the existence of polar BA L outflows (e.g. IZhou et all 
120061 IGhosh fcPunslvl 120081 ). This means that there is 
no simple orientation model which can explain all the 

features observed in BAL quasars. 

It has been suggested bv lBecker et a l. (2000) that, be- 
cause of their small sizes, most of the radio-loud BAL 
quasars belong to the class of compact radio sources, 
namely compact steep spectrum (CSS) objects and gi- 
gahertz peaked spectrum (GPS) objects. GPS and CSS 
sources are considered to be young radio sources with 
linear sizes less than 20kpc, entirely contained within 
the extent of the host galaxy. The high resolution ob- 
servation in the VLBI technique is the best way to 
learn about their morphologies and orientation. How- 
ever, very small fraction of compa ct BAL quasars have 
been observed with VLBI so far (|Jiang fc Wangl [20031 : 



iLiu et all 120081: 
About half of 
them still have unresolved radio structures even in the 
high resolution observations, the others have core-jet 
structures indicating some re-orientation or very complex 
morphology, suggesting a strong interaction with the sur- 
rounding medium. All of them are potentially smaller 
than their host galaxies. The analysis of the spectral 
shape, variability and polarization properties of some of 
them shows that they are indeed similar to CSS and GPS 
object s, and are not oriented along a particular line of 
sight ([Montenegro- Montes et al.ll2009D . 

1045+352 is a comp act radio-loud CSS source. The 
spectral observations (jWillott et al.l l2002t ) have shown 
1045+352 to be a quasar with a redshift of z — 
1.604. It has also been classified as a HiBALQSO 
based upon the observed very broad C IV absorption, 
and it is a very luminous submillimeter object with 
detec tions at both 850 /jm and 450 /xm (jWillott et al.l 
I2002T ). The radio luminosity of 1045+352 at 1.4 GHz 
is high (logLi.4GHz=28.25 W Hz -1 ), making this source 
one of the most radio-luminous BAL quasars, with 
a value similar to that of the first k nown radio- 
loud BAL QSO with an FRII structure (jGregg et al.1 
120001). The Chandra X-ray obse rvation of 1045+352 
( Kunert-Bairaszewska et all I2009D shows that the X- 
ray emission of 1045+352 is very weak in compari- 
son to the other radio-loud BAL quasars which, to- 
gether with the high value of the optical— X-ray index 
a ox — 1.88, suggests the presence of an X-ray absorber 
close/in the BLR region. The result of the spectral 
energy distribution (SED) modeling of 1045+352 indi- 
cate that the X-ray emission we observe from 1045+352 
may mostly be due to X-ray emission from the rela- 
tivistic jet, while the X-ray emission from the corona 
is absorbed in a large part. The first multifrcquency ra- 
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TABLE 1 
Basic parameters of 1045+352 



TABLE 2 

Flux densities of 1045+352 principal components from the 

5 GHz MERLIN and EVN images 



Parameter 



Value 



Source name (B1950) 1045+352 

RA (J2000) extracted from FIRST 10 h 48 m 34! 247 
Dec (J2000) extracted from FIRST +34° 57.'0 24'.'99 

Redshift z 1.604 

Total flux density Si.4GHz ( m Jy) 1051 

log.Li.4GH3 (W Hz" 1 ) 28.25 

Total flux density S4.85GHZ ( m Jy) 439 

Spectral index a\-f^^ 0.70 
Largest Angular Size measured 

in the 5 GHz MERLIN image 0'.'5 

Largest Linear Size (h" 1 kpc) 4.3 

Radio-loudness parameter R* 4.9(4.1) 

Intrinsic extinction Ay 1.5 

K-corrected absolute magnitude My -22.83(-24.33) 

Notes. Quantities in parentheses are corrected for intrinsic extinc- 
tion. Spectral index defined as (S oc i/~ Qr ) 

dio observations of the quasar 1045+35 2 were made by 
iKunert-Bairaszewska fe Mareckil ([2007D . The complex 
compact structure has been resolved into many subcom- 
ponents and indicates that the jet is moving in a non- 
uniform way in the central regions of the host galaxy. 
Here we present new, more sensitive, high resolution 
EVN+MERLIN 5 GHz observations of this source and 
analyze them. 

2. RADIO OBSERVATIONS 

1045+352 belongs to the primary sample of 60 can- 
didates for CSS sources selec ted from the VLA FIRST 
catalogue (jWhite et al.l I1997D . Initial observations of 
all the candidates w ere made with MERLIN at 5 GHz 
(jMarecki et al.l l2003t ) and 1.7, 5 and 8.4 GHz VLBA 
follow-up of 1045+352 was carried out on 13 Novem- 
ber 2004 in a snapshot mode with p hase-referencing 
(jKunert-Bairaszewska fc Mareckil l2007t ). The observa- 
tions showed that this source has a complex radio mor- 
phology with a probable radio jet axis reorientation. To 
look for the traces of the above-mentioned scenario we 
planned more sensitive high-resolution EVN+MERLIN 
5 GHz observation that was carried out on 2 June 2007 
in a full-track mode. The target source itself was ob- 
served for 7hrs. The whole data reduction process was 
carried out using standard AIPS procedures. IMAGR 
was used to produce the final images separately for EVN 
and MERLIN, and finally the combined EVN+MERLIN 
image was created (Fig. [1]). The flux densities of the 
main components of the target source were then mea- 
sured, by fitting Gaussian models, using task JMFIT 
(Table [2]). For more extended features the flux densi- 
ties were evaluated by means of IMSTAT (*, Table [2]). 
The basic parameters of 1045+452 have been gathered 
in TableHJ 

Throughout the paper, we assume the cosmology with 
H =71kms- 1 Mpc- 1 , fL M =0.27, ^a=0.73. 

3. RESULTS 

Our previous MERLIN image of 1045+352 
(|Kunert-Bairaszewska fc Mareckil [2007) shows this 
source to be extended in both NE/SW and NW/SE 
directions suggesting a reorientation of the jet axis. 

The new EVN+MERLIN 5 GHz radio observations of 
1045+352 BAL quasar revealed more details of its com- 
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plex morphology (Fig. [TJ. We identify the structure C 
with the radio core and the structures A - A 3 with the 
subsequent jet activity. The new 5 GHz EVN observa- 
tion shows a jet (indicated as a A, Fig.QJ,) emerging from 
the core (C) in the E/SE direction, and then another jet 
(Ai) emerging in the NE direction. The feature indicated 
as B is probably a counter-jet, and its trace was also 
visible in the previous lower resolution 1.7 G Hz VLBA 
image (jKunert-Bairaszewska fe Mareckil [20071. No cor- 
responding B1-B3 structures are found, possibly due to 
their weakness. There is an appreciable (~ 60°) mis- 
alignment between the axis of the inner structure (jet A) 
with respect to the outer structure (jet A\). The new 
combined image shows diffuse radio emission in the cen- 
ter of 1045+352 quasar suggesting very dense medium 
of the host galaxy of the source and indicating strong 
interactions (Fig. [Tb). 

The projected linear size of the whole struc- 
ture measured in our old snapsho t 5 GHz image 
(jKunert-Bairaszewska fe Mareckil |2007| ) as a size of the 
contour plot is 4.3 kpc (Table [J). The projected dis- 
tance between the radio core C and the feature A 3 
in the new 5 GHz MERLIN image (Fig. [TJd) mea- 
sured as a separation between the component peaks 
is 1.7 kpc. The projected size of the innermost part 
(A+C) is 190 pc. The size of the radio core amounts 
to 68 pc when measured in the 8.4 GHz VLBA image 
(jKunert-Bairaszewska fe Mareckill2007f ) where the inner- 
most components were resolved. The measured compo- 
nent fluxes from the 5 GHz EVN and MERLIN images 
are collected in Table [3] 

Based on the higher resolution VLBA image 
(jKunert-Bairaszewska fe Mareckil l2007h and using the 
radio -to-optical luminosity ratio ([Wills fe Brotherton I 
Il995h we have estimated the angle between the jet axis 
and the line of sight to be 9 ~ 30°. However, this should 
be treated as a rough a pproximation since the radi o-to- 
optical luminosity ratio ([Wills fc Brotherton 1119951 ) was 
defined for large scale objects and itself suffers from very 
large uncertainty. In addition, the relation we used is 
based on the radio luminosity of the core, so we estimate 
that the derived value of the angle corresponds to the 
innermost part of the jet - component A. 

Component A 3 (Fig. [TJd) seems to be the oldest struc- 
ture of 1045+352. A new emission - radio jet A is emerg- 
ing from the core C and an earlier phase of the jet activity 
is represented by components A\ and B. Because of their 
symmetry with respect to the radio core, we treated them 
as jet and counter-jet, respectively. The jet/counter-jet 
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MilliARC SEC 
Center at RA 10 48 34.2480005 DEC 34 57 25.041999 
Grey scale flux range- -2.0 223.3 mjy/beam 
peak flux density=223.29 mjy/beam, beam size=62.33 x 52.31 mas 
first contour level-0.35 mjy/beam 

Fig. 1.— Radio images of 1045+352: a) 5 GHz EVN image, b) 5 GHz MERLIN image, c) combined EVN+MERLIN 5 GHz image. 
Contours increase by a factor 2 and the first contour level corresponds to m 3<r. 



brightness ratio is given by: 



R 



1 + /3cos& 
1 - j3cos<d 



2+oy 



(1) 



where (3 is the jet/counter-jet speed, 9 is an angle be- 
tween the jet axis and the line of sight and cijet is the jet 
average spectral index. The value of the sideness param- 
eter is R ~ 13.5. Then assuming a typical value of the 
spectral index otj e t — 0.75 and the jet speed /? = 0.9 c 
we obtained the value of the angle between the jet axis 
and the line of sight to be 9 ~ 61°. However, 9 varies in 
the range 51° — 64° for j3 and ctj e t values in the range 
0.7c — 0.98c and 0.7 — 1.1 respectively. This may indi- 



cate a large change between the orientation of 1045+352 
radio structures on the different scales and fits well into 
the picture of complicated radio morphology. 

An assumption of 9 = 30° — 61° yields the deprojectcd 
linear size of the source of ~ 5 — 9 kpc, also indicating a 
young radio source. 

4. DISCUSSION 

The complex morphology and a change of the jet direc- 
tion observed in 1045+352 may result from (1) a merger, 
(2) a jet precession, or (3) jet-cloud interactions. Below, 
we discuss possible origins of the multiple radio struc- 
tures and their misalignment. 
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TABLE 3 
Models of the accretion disk evolution 
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Fig. 2. — The time evolution of the core luminosity for 2 values 
of mass and the mean accretion rate and 2 values of the viscosity 
parameter: m, = 0.36,0.68, M = 2 X 10 7 ,5.6 X 10 7 M Q , and a = 
0.03,0.1. The parameters are given in the Table[3] 



4.1. Quasar reactivation 
4.1.1. Galaxy merger 

Observations indicate, that about 50% of young AGN 
contain double nuclei in their host galaxies or exhibit 
morphological distortions that ar e supp osed to be due to 
the past merging events (|Q'Dealll99a iLiu 1 12004| ) . Al- 
though there are no obvious signs of merger event in the 
optical image of 1045+352 (from the Sloan Digital Sky 
Survey, SDSS), the radio distortions indicate this possi- 
bility. Therefore we discuss below the possible scenario 
of a merger event in QSO 1045+352. As shown by Liu 
(2004), the spin axis of the black hole formed after the 
merger changes its orientation from the vertical with re- 
spect to the outer accretion disk to the aligned with the 
rotation axis of the binary on timescale 10 5 yr. 

The angle by which the jet might have precessed in 
1045+352 is difficult to be determined observationally, 
due to orientation effects. The inclination angle of the 
A+C structure (jet axis) to the line of sight is about 30°. 
In the case of the older components A\ and B it seems 
to be rather 9 ~ 61°. The images of the radio structures 
are projected onto the plane of the sky. Therefore even 
for a very small precession angle, e.g. 10°, the misaligned 
structure will appear shifted by a large angle on the 2-D 
map, e.g. 90°, due to the projection effect. 

In the case of 1045+352, a timescale for the most dis- 
tant structure A3 (1.7 kpc, in projection) will be equal 
to t = 5.5 x 10 4 yrs if we assume the jet velocity of 0.1c, 
and t = 1.8 x 10 4 yrs for 0.3c. Therefore in this source, 
we may be witnessing an ongoing process of the disk re- 
alignment after a merger event, which took place after 
structure A3 was ignited. 

4.1.2. Accretion disk instability 
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Totes, a - viscosity parameter, T q - quiescent time, T ou t - outburst 



In this section, we apply the accretion disk instabil- 
ity model to explain the "reactivation" of the 1045+352 
core. Below, we describe in brief the basic properties of 
this model, which was described in mor e detail elsewhere 
(jJaniuk et alj|2002t ICzernv et al.ll2009ft . 

Instabilities of t he accretion disks have been shown to 
opera t e on long ( Siemiginowska et al. 1996; Janiu k et al.l 
l2004t Uaniuk k Czernvl I2007D and s hort timescales 
(| Janiuk et al.l 120021 : ICzernv et al.l 120091 ) , depending on 
their nature (the dwarf- nova type of instability, caused by 
the partial ionization of hydrogen; or the radiation pres- 
sure instability, studied in microquasars). Scaling the 
outbursts of the Galactic microquasar GRS 1915+105, 
lasting about 100—2000 seconds for a black hole mass of 
10M Q , to the quasar hosting a supermassive black hole, 



gives outbursts on timescales of 10 — 10 years. 

The unstable disk surrounding a black hole is subject 
to thermal and viscous instability if the radiation pres- 
sure dominates over the gas pressure. If the accretion 
rate outside the unstable region (i.e. the mean accretion 
rate) is such that the disk is in the unstable mode, the 
source enters a cycle of bright, hot states, separated by 
cold, quiescent states. In the hot state, the luminous 
quasar could power a radio jet, while during the cold 
state the radio activity ceases. 

The quantitative results, i.e. the outburst amplitudes 
and durations, are sensitive to the model parameters: 
black hole mass and viscosity coefficient. In addition, 
the description of the heating is essential. If the heating 
is proportional to the total pressure, the outburst am- 
plitudes are very large. The heating proportional to the 
square root of the gas times the total pressure reduces 
the amplitude of the disk outbursts. 

We compute the time dependent accretion disk model 
for several combinations of the main parameters: four 
values of the mean accretion rate, m = 0.36, 0.68, 0.7 and 
1.0 in units of the Eddington accretion rate; three values 
of the viscosity parameter: a = 0.01,0.03 and 0.1; and 
three values of the black hole mass, M — 2 x 10 7 , 5.6 x 10 7 
and 4 x 10 8 Af Q (see Table 3; models are labeled with 
A, B, and M). The accretion rate is in the units of the 
Eddington accretion rate, MEddi which corresponds to 
the Eddington luminosity. In Table [3] we give the re- 
sulting time for the quiescent and outburst periods. In 
Figure [2] we show the exemplary light curves from the 
time dependent accretion disk model. In general, the 
ou tburst separation is consisten t with the one required 
by iRevnolds fc Begelman (1997), who proposed a phe- 
nomenological model to fit the observed number of radio 
sources and their sizes. 
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In our model, for a given black hole mass, the larger 
the mean accretion rate, the longer the duration of a cy- 
cle episode, both in hot and cold states. In addition the 
a parameter affects the results and for smaller viscosity, 
th e cycle duration is lon ger. This relation was calibrated 
by ICzernv et al.l (|2009l ) using a grid of models for vari- 
ous black hole masses and Eddington ratios. The proper 
value of a bolometric correction in this relation can give 
us an estimation of an upper limit for the age of a radio 
source (i.e. the duration of the hot state) . The predicted 
timescales and durations are in agreemen t with observa- 
tions of compact radio sources (|Wull2009T ). 

The outbursts are associated with the ejections of ra- 
dio jets. The jets are then turned-off between the out- 
bursts and each radio structure will represent a new out- 
burst. In case of an apparently young, compact source 
we can suspect that in fact it is an old, reactivated ob- 
ject, in which the vast radio structures have already faded 
away and are not visible. This mechanism would explain 
the apparent statistical excess of the compact sources 
with r espect to the ga l axies with extended radio struc- 
tures (jO'Dea fc Bauml (|1997l )). On the other hand, the 
timescales for fading of a radio source may be much 
longer than the separation timescales between the out- 
bursts and we should observe the fossil radio structures in 
addition to the compact one (|Baum et al.lll990l) . These 
kind of morphologies have already been observed in large 
scale ob jects showing evidence of two or more active per i- 
ods (e.g lSchoenmakers et al.ll2000l : Uamrozv et ai1l2007j ). 

Here, we present the results of the modeling of the 
intermittent activity for the source 1045+352. The 
best fitting parameters for the q uasar 1045+352 from 
iKunert-Bairaszewska et all ()2009l ) were M BH = 2 x 
1O 7 M , to = 0.36 and a = 0.1 (model Ml). Such a large 
accretion rate is high enough for the radiation pressure 
instability to operate (M cr it ~ 0.2MEdd)- We take these 
parameters as an input to our time-evolution computa- 
tions, and for comparison we calculate model M2, with a 
different viscosity. The accretion rate in solar masses per 
year is therefore M = to x 3.52 M 8 = 0.253 M 8 , where M$ 
is the black hole mass in the units of 10 8 Solar masses. 
Because the black hole mass determination for this source 
is not very certain, and seems to be rather a lower limit, 
we calculated also several models for somewhat larger 
masses (models M3, M4, for which we used the mass es- 
timation from the monochromatic luminosity; models A 
and B, where the assumed mass value enables the radio 
source to escape from the host galaxy). The model is 
quite sensitive for the adopted central mass value and 
the results presented here are only examples. Full pa- 
rameter space for the mass and ac cretion rates for larg e 
sample of quasars was presented in lCzernv et al.l (J2009). 

In the bottom left panel of Fig. [2] we plot the 
light curve for the given best fit parameters, taken 
from the spectral modeling ()Kunert-Bairaszewska et al.l 
2009). The parameters are given in the Table [3J model 
Ml. The resulting duration of an active phase in this case 
would be about 25 years , while for the quie scent phase 
it is about 45 years (see ICzernv et al.l l2009f ) for the de- 
pendence between the activity timescale and bolometric 
luminosity of the source) . 

In this scenario, the subsequent radio structures would 
be fed by the jet produced during the active phase of 



the central engine. In the quiescent phase, the jet stops 
feeding the radio lobe, and further expansion of the rem- 
nant radio structure proceeds due to the accumulated jet 
energy. If the proper motion of the outer jet components 
in 1045+352 is in the order of 0.1 — 0.2 c, and the restart 
of the activity was after T q ~ 50 yrs, then the physi- 
cal separation of the radio structures would be of about 
10 pc. If the velocities of the jet were rather of 0.7 — 0.9 c, 
such as those characteristic of the inner structures A\ - 
B, then this separation would increase to 35— 45 pc. 

During the active phase, the radio lobes are fed by the 
expanding jet and the source grows. For the outburst 
duration calculated in model Ml, the size of the central 
radio core would be only about 1.5 pc for a jet velocity of 
0.1 c, 4.6 pc for velocity of 0.3 c and 15 pc for a relativistic 
jet speed of 0.9 c. The resulting size of the radio source 
is therefore too small in comparison with the structure 
size estimated from our images (the size of the radio core 
is ^70 pc). The much larger actual size requires longer 
activity timescales, which are possible in the models with 
larger black hole mass and accretion rate and/or smaller 
viscosity. We therefore checked a number of additional 
models, with larger black hole mass and accretion rate. 
In particular, the subsequent models M2, M3 and M4 
adopt these parameters within the range of the plausible 
values for our previous spectral modeling. These values 
were not the best fit, however, they also gave a reasonable 
shape of the high energy spectrum of this quasar and 
therefore cannot be excluded. As Table [3] shows, model 
M4, with black hole mass of 5.6 x 10 7 M Q , accretion rate 
of to = 0.68 and viscosity a — 0.03 results in an outburst 
timescale of almost 900 years, which gives a size of the 
radio core of 30—80 pc, depending on the jet velocity 
(non-relativistic) . This is in good agreement with the 
estimated size of the central radio core. On the other 
hand, if the jet velocity is relativistic, the models with 
smaller black hole mass would be favored, consistent with 
the spectral fitting results. 

The above model of the cyclic activity of the central 
core is based on the 1-D, vertically averaged calculations 
and does not describe any effects in the polar direction. 
Therefore the jet axis assumed to be perpendicular to 
the disk plane is constant in time, and the subsequent 
radio structures should be co-aligned. Taking into ac- 
count, that the lobe separation is much smaller than the 
distance to the source, the viewing angle of 9 ~ 30° be- 
tween the line of sight and the younger radio lobe would 
be practically the same as the angle to the older lobe. 
In the case of 1045+352 the parts of the jet we observed 
in the 5 GHz structure are not co-aligned. Structures 
A\ - A2 are shifted with respect to the core, while the 
structure A3 is misaligned by a large angle. 

4.2. Precession of the jet 

The observed misalignment between the young and old 
radio structures could be the result of the changed direc- 
tion of the jet axis between the activity episodes. The 
precession of the jet could be caused by various mecha- 
nisms, such as tidal forces, disk warping, irradiation , and 
the Bardeen-Petterson effect (|Caproni et al.ll2006|) . In 
this section, we consider the precession of the innermost 
accretion disk due to interna l instabilities within the ac - 
cretion flow. As studied by Uaniuk et all ()2008l , 12009ft . 
acoustic instabilities of the Papaloizou-Pringle type can 



G 



M. Kunert-Bajraszewska et al. 



arise in the highly supersonic inner parts of the accretion 
flow. The azimuthal modes of such instabilities will lead 
to the initially differential tilt of the rotating torus and 
its subsequent precession. The 3-D hydrodynamic sim- 
ulations were performed to determine the rate of such 
precession, which is dependent on the adopted equation 
of state (i.e. the gas adiabatic index), and scales with 
the black hole mass. 

In Figure [3] we show the precession angle as a function 
of time, for Mbh = 2 x 10 7 M©. The simulations were 
performed for a few equations of state, and we show here 
only the result for the adiabatic index 7 = 5/3, which 
corresponds to the gas pressure dominated case. The 
twist angle, shown in the Figure as a function of time, is 
defined as follows: 
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arccos 



\ T1 1 Til 
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where L x , L y and L z are the Cartesian components of the 
angular momentum vector. The twist is therefore defined 
as a cumulative angle by which the angular momentum 
vector revolves in the x — y plane by the time t. 

Initially, when the torus is axisymmetric and no tilt 
occurs, the only component of the angular momentum 
vector is L z and its sign depends on the direction of 
the flow rotation. When the torus tilts, the non-zero 
L x and L y components appear and the rotation axis tilts 
toward the x — y plane. Initially, the torus precession is 
differential and the twist angle rises fast mostly in the 
innermost parts of the flow. In the Figure [3l we show 
the twist averaged over radius between the inner edge 
and 50 Schwarzschild radii. The twist angle decreasing 
from ~ 60° to ~ —80° in the innermost torus reflects 
the clockwise precession. The moment, when the preces- 
sion starts during the hydrodynamic simulation, as well 
as the period of precession, depend on the model (i.e. 
on the equation of state). For the model shown here, 
we estimated that the precession period is 2.4 x 10 4 in 
the units of dynamical time at the inner radius. For the 
black hole mass of 2 x 10 7 M Q it is about 2.5 years 

The simulations we described above assume a radia- 
tively inefficient flow, and therefore they are adequate 
for the quiescent state of the central core of the quasar. 
In other words, the accretion rate in the processing disk is 
negligibly small, and as such is not taken into account in 
the 3-D modeling. Accordingly, we postulate here that 
the precession of the disk occurs between the activity 
episodes, and each subsequent episode restarts the jet 
production along a misaligned axis. The direction of the 
axis is basically random, because the period of preces- 
sion is much smaller (about 20 times) than the duration 
of the quiescent phase. In consequence, the younger radio 
structure is easily formed at a different angle to the line 
of sight than the older one and the probability of that is 
very high taking into account the very small precession 
period. 

We note that this approach is a simplification and in or- 
der to really model the precession effect we would need to 
compute the 3-dimensional model with somewhat higher 
accretion rate and radiative efficiency. This is beyond 
the scope of our present co mputations. However, based 
on the results of lDasI ()2004f ). we suppose that at least for 
not extremely large accretion rates and viscosities, the 
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Fig. 3. — The twist angle averaged over radius, from 1.5 to 50 
.Rschwi as a function of time. The black hole mass is 2 X 10 7 Mq. 

effect of precession would be enhanced, because of the 
outward shifting of the sonic point. 

4.3. Interaction with the interstellar medium, 

The repetitive outbursts of the central core due to the 
accretion disk instability may potentially have a strong 
impact on the ISM. Due to the cyclic irradiation from the 
central UV source, the shocks in the ISM may arise on 
short timescales and keep the gas warm. The hot spots 
can travel to a certain distance within the host galaxy, 
before the outburst is finished. This depends on the out- 
burst duration and in case of a short timescale the radio 
source doe s not expand beyon d the host galaxy. Ac- 
cording to ICzernv et al.1 (J2009T ) in order for the source 
to escape the host galaxy the outburst would need to 
be longer than 10 4 yrs. For 1045+352, however, this 
would require a large black hole mass in this quasar, 
M = 4 x 1O 8 M , an accretion rate equal to the Ed- 
dington limit, and a very small viscosity (models A and 
B in Fig. [3]). Such values were ex cluded by our prev ious 
spectral modeling (jKunert-Bajraszewska et aLll2009l ). In 
particular, for 1045+352, the estimated core size is about 
70 pc. 

After the jet activity in 1045+352 ceases the radio 
structure will be driven further due to the pressure ex- 
pansion, up to a distance of about I kpc (for the ISM 
density of n = p/m p — 0.1 cm -3 , its temperature 
of 10 7 K and the jet power of L, Rt = 10 44 erg s _1 
(|Kunert-Bairaszewska fc Ma recki 2003)), assuming that 
the accumulated jet energy is comparable to the en- 
ergy content in the h eated medium, E = R 3 pc^ ~ L ict t 
(jStawarz et al.l f2008h . This distance is comparable to 
that of the outermost radio structure, A3 in our map 
(1.7 kpc). If the jet power was somewhat smaller, the 
distance of the radio structures driven by the pressure 
expansion would also be smaller, and could be compara- 
ble to the inner component Ai as well. 
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When the expansion stops, the radio structure will 
start recollapsing. The recollapse time is typically much 
longer than the outburst cycle, and the "tunnel" made 
by the expanding jet will not close before the next out- 
burst. (For our favorable parameters, the next outburst 
will start after about 4 x 10 3 years.) Therefore the jets 
will always propagate into a rarefied medium. This would 
imply that the jet is propagating with a large veloc- 
ity, and thus t he velocity of > 0.1c would be favored 
(Stawarz 2004). In this case, any distortions in the jet 
morphology would propagate rather quickly and the ap- 
parent twist seen in the maps (i.e. structure A3) could 
be enhanced. It has to be noted however, that there is 
a large uncertainty in the measurements of a black hole 
mass. According to the SED modeling the value of the 
black hole mass in 1045+352 is on the order of 10 7 M Q , 
according to the Mgll A2800FWHM it is on the order 
of 10 8 M Q . We therefore cannot exclude the possibility 
that the radio structure of 1045+352 will be able to es- 
cape from the host galaxy during this phase of the radio 
jet activity. 

Our analysis of radio components observed in 
1045+352 may suggest the following interpretation: the 
radio jet emerging from the core in a E/SE direction is 
not able to get through the dense environment and bends 
to NE direction (component A\). The radio emission in- 
dicated with Ai (Fig. [TJ;) and visible in a NE direction 
is probably another trace of jet-cloud interaction. Its 
continuation (A3) is visible on the MERLIN 5 GHz im- 
age (Fig. |Tp). The weakness of the counter-jet emission 
is probably caused by the large beaming. However, the 
structure may be as well interpreted as involving at least 
three phases of quasar activity: components A2-A3 as 
the oldest one, structure Ai-B as the younger one, and 
the jet A as the current activity direction. 

4.4. Radio structures of BAL quasars 

Observations based on the VLA FIRST survey have 
revealed the existence of radio-lou d BAL quasars 
(jBecker et aJj|2000t iMenou et all 120011 ), which together 
with th e 5 radio-loud BAL quas ars from NVSS discov- 
ered by iBrotherton et al.l (|1998f ) make a large popula- 
tion of objec t s, ab out 14%-18% of the total sample of 
iBecker et al.1 (|2000l ). There are only ten known radio- 
loud BAL qu asars with the extende d doublc-lob ed ra- 
dio emission (IWills et aUll999t iGregg et al.1 [2000. 2006; 
IBrotherton et al.ll2002t IZhou et al.ll2006D . Among these 
are: FRII quasar with a ste ep spectrum core s uggest- 
ing restarted activity inside (jGregg et al.ll2000l) . a hy- 
brid object whose morphology indicates some type of 
jet-medium interactions (jWills et al.l I1999T ). and a very 
core-dominated radio t riple, probably beamed source 
(|Brotherton et al.l l2002f ). However, most of the BAL 
quasa rs remain unresolved even with the VLBI observa - 
tions (jMontenegro-Montes et all 120091 : iDoi et al.ll2009f ). 
This implies they are either beamed with the jet point- 
ing close to the line of sight, or they are very compact 
objects starting their activity. The resolved ones show 
complex double or triple, in many cases asymmetric, 
morphology. In general, the radio-loud BALs tend to be 
compact in the radio, similar to GPS and CSS sources, 
which are thought to be the y oung counterparts of pow- 
erful large-scale radio sources (|Becker et al.ll2000l) . They 
have wide range of spectral indices which together with 



the polarization properties suggest that the radio-loud 
BAL qua sars are not oriented along a particular line 
of sig ht ([Becker et al.l 120001: iMontenegro-Montes et al.l 
2009), contrary to the orientation model. The evolution 
interpretation of the nature of BAL phenomena is also 
still discussed. 

1045+352 is a radio-loud BAL quasar and CSS object 
with a very complex structure which suggests that the 
radio jet is oriented close to the line of sight. Morpholo- 
gies like th at observed in 1045+352 are seen a mong other 
CSSs (e.g. iFanti et~afl 12001 lAn et al.l l2010l) . although 
they are not very common. It is very likely that many 
CSS objects interact with an asymmetric medium in the 
central regions of their host galaxies, and t his can cause 
the observed as ymmet ries an d distortions (Saiki a et al.l 
12001b Uevakumar et al.ll2005[ ). In the case of 1045+352 
an additional fact may matter, namely, thar our Chandra 
X-ray ob servations suggest a presence of t he absorbing 
material (jKunert-Bajraszewska et aLll2009l ) close/in the 
BLR region. It is then possible that the outflowing mate- 
rial detected in absorption may also have interacted with 
the radio jet disturbing its structure (see also the inter- 
pretation of the rad io morphology of CSS quasar 3C48 by 
iGupta et al.l (|2005[) ). Such a strong interaction between 
a relativi stic jet and a BAL w ind has been reported in 
Mrk231 ([Reynolds et al.ll2009t ). although the interaction 
takes place at the small deprojected distance of ~ 4 pc 
from the r adio core. 

Recently IShankar at al.1 (J2008I ) studied the radio prop- 
erties of BAL quasars based on the SDSS/DR3 and 
FIRST surveys and found that the number of classical 
BAL QSOs decreases with increasing radio power. They 
discussed several plausible physical models which may 
explain the observations, namely the evolution and ori- 
entation model, however none of them can fully explain 
all the features and correlations observed among radio- 
loud BAL quasars. As has already been suggested (jElvisI 
2000) not only equatorial disk wind but also cone out- 
flows can be responsible for BAL features. Thus an in- 
ternal orientation effect such as the opening angle of the 
accretion disk wind instead of the external orientation ef- 
fect such as an angle between the jet axis and our line of 
sight, could be what determines a possibility of observing 
any BAL features. 

The discussed orientation of 1045+352 is not in 
agreement with the orientation model, however, the 
BAL features are clearly visible in this source. After 
iKunert-Bairaszewska et all (J20091 ) we suggest that there 
is no a direct dependence between the radio jet orienta- 
tion and the probability of BAL visibility. 

5. SUMMARY 

1045+352 is a CSS object and a HiBAL quasar at a 
medium redshift. The spectrum classification and linear 
size of the source indicate that it is a young object in 
the early phase of its evolution. The new more sensi- 
tive high-resolution EVN+MERLIN 5 GHz observations 
revealed details of the source radio emission and con- 
firmed its complicated radio structure. The radio mor- 
phology of 1045+352 is dominated by the strong radio 
jet resolved into many subcomponents and changing its 
orientation during propagation in the central regions of 
the host galaxy. The presence of the dense environment 
and material outflows have been already confirmed in 
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this source based on the submillimeter and X-ray obser- 
vations and as we discussed, they can have a significant 
impact on the radio jet morphology of 1045+352. We 
conclude that the most probable interpretation of the ob- 
served radio structure of 1045+352 is the ongoing process 
of the jet precession due to internal instabilities within 
the flow and/or jet interaction with a dense, inhomoge- 
neous interstellar medium. 

It is difficult to establish the orientation between the 
jet axis and the observer in 1045+352 because of the 
complex structure. Nevertheless taking into account the 
innermost radio structure we suggest that the radio jet 
is oriented close to the line of sight. This means that the 
opening angle of the accretion disk wind can be large in 



this source. 
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